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TTic  study  reported  herein  was  arnducted  within  the  facilities  of  the  Envi¬ 
ronmental  Laboratory  (EL)  of  the  L'.S.  Army  Engineer  Waterways  Experiment 
Station  :  \VES)  during  the  period  of  May  1991  to  June  1992  for  Headquarters, 
U.S.  Army  Corps  of  Engineers. 
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Harrison,  Director,  EL. 
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and  prepared  the  report  under  the  direct  supervision  of  Messrs.  Malcolm 
Kcown  and  Kertneth  Hall,  Chiefs.  Environmental  Constraints  Group,  ESD. 

At  the  time  of  publication  of  this  report.  Director  of  W'ES  was 
Dr.  Roben  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 

This  report  should  be  cited  as  follows: 

Curtis.  J.  O.  (1993).  "Microwave  dielectric  behavior  of  soils; 

Report  2,  A  unique  coaxial  measurement  apparatus,”  Technical 
Repon  EL-93-25,  U.S.  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 
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Introduction 


New  and  improved  methods  of  remote  sensing  have  increased  the  under¬ 
standing  of  Earth's  origins,  its  resources,  and  those  processes  that  contribute 
to  its  dynamic  (on  a  large  time  scale)  nature.  Scientists  and  engineers  from 
many  disciplines  are  constantly  exploring  new  ways  to  quantify  the  Earth's 
propenies  for  their  particular  applications.  Among  these  methods  are  mea¬ 
surements  of  electromagnetic  energy  in  many  different  wavelength  regimes, 
both  passive  and  active.  The  development  of  small  powerful  sources  and 
ultrasensitive  receivers  along  with  improved  data  processing  capabilities  has 
fostered  renewed  interest  in  measurements  of  the  Earth  within  the  microwave 
region  of  the  spectrum  where  wavelengths  in  air  range  from  a  few  millimeters 
to  several  meters, 

Namral  terrain  surfaces  consist  of  bare  soils,  rocks,  vegetation,  and  water. 
This  study  was  initiated  with  the  hope  of  making  a  meaningful  contribution  to 
the  understanding  of  microwave  interactions  with  natural  terrain.  The  first 
report  in  this  series  attempts  to  summarize  what  is  currently  known  about  the 
measurement  and  modeling  of  the  electrical  properties  of  well-characterized 
soils  (Curtis,  In  Preparation).  This  report  contains  a  description  of  a  new 
dielectric  property  measurement  capability  recently  developed  at  the 
U.S.  Army  Engineer  Waterways  Experiment  Station  in  Vicksburg,  MS,  that 
will  be  used  in  the  future  to  establish  a  comprehensive  database  of  soil  proper¬ 
ties  covering  a  broad  spectrum  of  soil  types,  measurement  frequencies,  sample 
moisture  contents,  and  sample  temperatures.  The  final  report  of  the  series 
will  present  some  initial  measurement  results  and  will  discuss  two  modeling 
approaches  to  interpreting  those  data.  Most  of  the  information  contained  in 
this  series  of  reports  is  drawn  from  a  doctoral  research  program  recently  com¬ 
pleted  by  the  author  (Curtis  1992). 

In  an  effort  to  make  each  report  of  this  series  as  self-contained  as  possible, 
the  rationale  for  pursuing  this  research  program  as  well  as  some  supporting 
information  in  appendices  is  included  within  each  document.  It  is  hoped  that 
this  will  be  useful  for  the  individual  who  might  read  only  one  of  the  reports. 
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Reasons  for  Studying  Electrical  Properties  of  Soils 


Soil  moisture  measurements 

The  dominant  factor  that  controls  the  electrical  behavioi  of  soils  is  the 
presence  of  water  (Topp,  Davis,  -md  Annan  1980).  Obviously,  one  would 
hope  to  take  advantage  of  this  experimental  fact  to  develop  a  means  of  accu¬ 
rately  and  quickly  measuring  the  moisture  in  soils  without  having  to  collect 
numerous  samples  in  the  field,  weigh  them,  dry  them  for  extended  periods  of 
time,  and  weigh  them  again  to  obtain  either  a  gravimetric  (weight  of 
w'ater/weight  of  dry  soil)  moisture  content  or  a  more  useful  volumetric  (vol¬ 
ume  of  water/ volume  of  soil  sample)  moisture  content. 

Numerous  attempts  have  been  made  to  develop  a  useful  method  for  mea¬ 
suring  soil  moisture  content.  All  have  met  with  varying  degrees  of  success, 
and  none  has  proven  accurate  under  all  conditions.  For  example,  a  technique 
that  essentially  amounts  to  bury  ing  radar  transmit  and  receive  antennas  in  the 
soil  and  relating  the  measurements  of  attenuated  received  signals  to  moisture 
content  (Birchak  et  al.  iQT'J)  is  destructive  to  the  soil  fabric  (the  way  in  which 
soil  particles  are  arranged),  is  very  much  controlled  by  the  fabric  and  the  size 
distribution  of  panicles,  and  precludes  the  use  of  the  same  instrument  in  mul¬ 
tiple  locations  (nature  is  not  homogeneous)  as  well  as  the  ability  to  easily 
repair  defective  equipment. 

Another  approach  taken  by  some  researchers  for  making  field  measure¬ 
ments  of  soil  moisture  (that  also  has  numerous  application  in  the  biomedical 
field)  is  that  of  measuring  the  change  in  fringe  capacitance  of  an  open-ended 
coaxial  probe  (Thomas  1966;  Brunfeldt  1987;  Gabriel,  Grant,  and  Young 
1986).  When  pressed  against  a  soil  whose  propenies  are  unknown,  the  result¬ 
ing  change  in  capacitance  produced  by  the  impedance  mismatch  is  related  to 
electrical  properties  through  calibration  relationships.  Problems  arising  from 
these  measurements  include  the  need  to  have  proper  contact  between  the  probe 
tip  and  the  soil  surface,  the  fact  that  the  volume  of  material  associated  with 
the  fringe  capacitance  is  quite  small  (1  cm’  or  less),  and  that  calibration  condi¬ 
tions  simply  cannot  account  for  all  of  the  dielectric  loss  mechanisms  that  exist 
in  natural  soils.  The  losses  in  moist  soils  can  be  highly  frequency  dependent 
over  a  range  of  several  frequency  decades  on  the  electromagnetic  spectrum. 

A  recent  variation  on  the  open-ended  probe  measurement  scheme  involves  the 
use  of  a  waveguide  section  instead  of  a  coaxial  device  (Parchomchuk, 
Wallender,  and  King  1990). 

The  concern  over  small  sample  volumes  can  be  overcome  with  a  redesign 
of  the  open-ended  coaxial  probe  that  replaces  the  solid  outer  conductor  with 
several  pointed  tines  (Campbell  1988)  which  allows  for  the  probe  to  be  pushed 
into  the  surface  of  soft  soils.  The  volume  of  soil  enclosed  by  such  probes  can 
easily  be  tens  of  cubic  centimeters.  The  tined  coaxial  probes  that  have  been 
built  to  date  operate  in  a  frequency  range  that  is  very  much  subject  to  the 
material-dependent  loss  mechanisms  alluded  to  above.  Research  continues  that 
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is  directed  to  better  understanding  these  various  mechanisms  and  to  fabricate  a 
probe  that  operates  in  a  frequency  range  that  is  not  subject  to  sucli  material- 
dependent  anomalies. 

Subterranean  investigations 

There  are  many  t>pes  of  electrical  measurements  m  soils  (and  rocks)  that 
are  used  to  understand  what  lies  beneath  the  Earth's  surface  (Telford  et  al 
1984).  Among  these  are  various  schemes  for  measuring  soil  resistivity  as  well 
as  the  attenuation  of  propagating  electromagnetic  waves 

Resistivity  measurements.  Resistivity  data  in  soils  are  collected  by  inject¬ 
ing  known  currents  (usually  of  a  frequency  less  than  60  Hz.  alternating  cur¬ 
rent  required  to  minimize  effects  of  charge  buildup  on  the  probe)  into  the 
ground  and  measuring  potential  differences  across  pairs  of  nearby  electrodes. 
Assuming  homogeneous  media  and  uniform  resistivity,  it  is  possible  to  calcu¬ 
late  from  the  potential  differences  an  apparent  resistivity  of  the  earth  material. 
These  measurements  will  be  affected  by  the  presence  of  water  in  the  soil  or 
rock,  by  the  presence  of  mineral  compounds  that  could  go  into  ionic  solution 
with  available  water,  and  by  the  physical  structure  of  the  subsurface  terrain 
itself  This  being  the  case,  resistivity  measurements  are  useful  as  a  measure 
of  subsurface  water  volume,  the  locations  of  mineral  deposits,  and  subsurface 
structure. 

Electromagnetic  wave  propagation.  Another  method  of  making  electrical 
subsurface  measurements  involves  the  transmission  of  electromagnetic  waves 
into  the  soil  and  detection  of  energy  that  results  from  waves  reflected  from 
subsurface  anomalies.  Having  the  ability  to  track  wave  propagation  in  time, 
either  by  pulsing  the  source  or  sweeping  over  a  known  frequency  band  in 
some  controlled  manner,  means  that  electromagnetic  wave  propagation  meth¬ 
ods  of  subterranean  investigations  are  particularly  useful  for  locating  the  depth 
of  electrical  anomalies  such  as  the  water  table  in  sandy  soil  (Olhoeft  1983; 
Stewart  1982;  Wright,  Olhoeft.  and  Watts  1984),  buried  pipes  or  wires,  or 
cavities  such  as  tunnels  or  caves  (Ballard  1983).  Other  applications  include 
the  delineation  of  stratified  media  (Lundien  1972'  and  determination  of  the 
thickness  of  ice  and  frost  layers  (Jakkula,  Ylinen,  and  Tiuri  19c O'Neill  and 
Arcone  1991). 

There  are  some  practical  bounds  on  the  utility  of  radio  frequency  systems 
to  conduct  subterranean  investigations  because  of  the  phenomenon  of  "skin 
depth,"  a  measure  of  the  attenuation  of  the  electromagnetic  energy  as  it  travels 
through  the  medium.  For  relatively  low-frequency  sources  (approximately 
200  MHz),  it  can  be  shown  that  low-loss  soils  such  as  dry  sands  can  possess  a 
skin  depth  of  approximately  10  to  15  m,  while  high-loss  soils  such  as  wet  silts 
and  clays  may  have  skin  depths  of  only  a  few  centimeters.  Modern  radio  fre¬ 
quency  receivers  arc  extremely  sensitive  devices,  often  having  a  dynamic 
range  of  50  to  100  db.  The  signal  at  skin  depth  represents  about  an  8.7-db 
loss  in  power  or  a  two-way  loss  at  the  receiver  for  reflected  signals  of  about 
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17.4  dh.  It  is  ccnainly  not  inconceivable  that  radio  frequency  receivers 
sliould  be  capable  of  successfully  detectinc  reflected  signals  from  subterranean 
anomalies  at  depths  of  two  to  three  skin  depths  or  more. 


Remote  sensing  of  environment 

Vinually  all  remote  sensing  of  the  environment  from  airborne  or 
spacebome  platfomis  involves  the  measurement  of  electrcmagnetic  radiation 
from  the  Eaith's  surface  and'or  atmosphere.  Passive  surveillance  involves 
measurements  of  emitted  radiation  and  that  reflected  from  natural  sources  such 
as  the  sun,  the  atmosphere,  and  surrounding  terrain  Active  remote  sensing 
measurement  systems  include  a  source  for  illuminating  the  target  of  interest. 
Whether  passive  or  active,  visual,  thermal  infrared,  microwave,  or  millimeter 
wave  sensors  are  utilized,  remote  seasing  is  the  collection  and  interpretation  of 
electromagnetic  radiation  and,  as  such,  demands  an  understanding  of  the 
dielectric  properties  of  those  materials  being  observed 

Environmental  remote  sensing  a''plications  form  a  list  that  grows  yearly  as 
electronic  components  arc  improved  and  data  collection  and  processing  hard¬ 
ware  and  software  become  faster,  more  reliable,  and  less  cosily.  For  exam¬ 
ple,  satellites  can  provide  worldwide  surveys  of  land-use  paiicms  to  monitor 
the  threat  of  urbanization,  waste  disposal,  and  erosion  of  the  land  (Colwell 
1983).  Similar  systems  (including  those  mounted  in  aircraft)  can  monitor  the 
health  of  vegetation  to  keep  abreast  of  such  things  as  lo.ss  of  forest  and  the 
potential  for  food  shonages.  Sea  traffic  in  the  far  northern  and  southern 
shipping  lanes  can  be  made  safer  through  the  use  of  airborne  and  spacebome 
Sensors  to  detect  ice  hazards. 

One  of  the  more  obvious  applications  for  microwave  remote  sensing 
devices  is  that  of  conducting  surface  moisture  surveys  to  help  predict  ground- 
water  availability  and  the  potential  for  flooding  Attempts  have  been  made  to 
relate  soil  moisture  to  both  laboratory  leflectance  data  (Lundien  1966)  and  to 
airborne  sensor  backscatter  measurements  (Ulaby,  Cihlar,  and  Moore  1974, 
John  1992)  Careful  airborne  sensor  measurements  might  provide  a  first 
approximation  to  the  complex  dielecf  .c  constant  of  the  soil  near  the  surface. 

Because  of  the  difference  in  dielectric  behavior  of  liquid  water  and  various 
forms  of  ice,  it  may  be  possible  to  use  airborne  sensors  to  detect  freezing  and 
thawing  (Wegmuller  1990)  in  remote  locations  that  could  be  used  to  predict 
spring  runoff  conditions  and  all  that  that  encompasses  for  agricultural  applica¬ 
tions,  the  effects  on  the  fishing  industry,  and  the  anticipation  of  flooding  in 
built-up  areas.  Military  analysts  are  concerned  about  soil  moisture  conditions 
because  of  its  impact  on  trafficability.  the  ability  of  vehicles  to  move 
effectively  over  natural  terrain. 

Another  remote  sensing  application  is  the  mapping  of  exposed  soils  and 
rocks  in  remote  areas  of  the  world  from  high-flying  aircraft  or  satellites, 
which  might  prove  useful  for  geomorphological  studies  (Swanson  et  al.  1988) 
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or  even  mineral  cxplora'.'on.  A  less  obvious,  but  recent,  application  of  micro- 
wave  remote  sensing  in  soils  dealt  with  archeological  surveys  in  desert  areas 
(Berlin  et  al.  1986;  McCauley  ct  al  1986). 


Others 

While  the  above  paragraphs  emphasize  some  of  the  most  obvious  and  use¬ 
ful  applications  of  a  better  understanding  of  soil  electrical  properties,  others 
have  been  noted  in  the  literature.  For  example,  some  researchers  have 
attempted  to  relate  electrical  property  measurements  to  the  physical  properties 
of  soils  (Campbell  and  Ulrichs  1969;  Hayre  1970,  .Arulanandan  and  Smith 
1973;  Madden  1974).  Of  course,  nothing  has  been  said  about  the  military’s 
need  to  better  understand  the  microwave  response  of  soils  that  form  the  back¬ 
grounds  to  militar)  targets,  i.e.,  when  and  why  does  clutter  become  a  source 
of  target-like  signatures"’  Outside  of  the  topic  of  soils,  studies  of  the  micro¬ 
wave  respon.se  of  foodstuffs  has  direct  application  to  quality  conirol  concerns 
in  the  food  industry  (Nelson  1973.  Nelson  1983). 

Another  new  application  of  this  technology  that  is  closely  related  to  the 
discussion  on  soil  moisture  is  that  of  detecting  liquid  ground  contaminants, 
either  near  the  surface  or  at  arbitrary  depths  using  a  specially  fabricated 
probe.  If,  as  wall  be  argued  later,  polarizable  liquids  can  be  characterized  by 
a  unique  frequency  of  peak  losses  because  of  the  dielectric  relaxation  phenom¬ 
enon,  then  a  probe  could  be  designed  to  measure  losses  over  a  frequency  span 
that  IS  broad  enough  to  detect  a  peak  loss  frequency  and,  coupled  with  the 
results  of  a  thorough  experimental  program,  to  identify  the  particular  contami¬ 
nant. 


Problem  Statement 


The  overall  objective  of  this  research  is  to  measure  and  model  the  dielec¬ 
tric  response  of  moist  soils  in  w'ays  that  will  test  the  current  understanding  of 
loss  mechanisms  over  a  broad  range  of  moisture  levels,  frequencies,  and 
material  temperatures  and  perhaps  suggest  new  ways  of  thinking  about  losses 
as  a  function  of  these  variables. 

■What  is  first  required  to  help  achieve  this  objective  is  an  apparams  that  will 
allow  for  the  measurement  of  attenuation  in  moist  soils  over  a  broad  spectrum 
of  frequencies  w'hile  maintaining  a  required  sample  temperature.  A  methodol¬ 
ogy  must  exist  for  allowing  sample  moisture  to  vary  in  a  measurable  way. 

The  remainder  of  this  report  describes  such  a  device  and  a  methodology  and 
includes  measured  data  that  substantiates  the  algorithm  used  to  compute  values 
of  the  complex  dielectric  constant  for  each  set  of  variables 
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2  A  Coaxial  Apparatus  for 
Dielectric  Measurements 


The  choices  for  laboratory  techniques  for  making  material  electrical 
propeny  measurements  include  resonant  cavity,  waveguide,  and  coaxial,  with 
some  new  work  being  done  on  hybrids  of  the  latter  two  (Taherian  ct  al. 

1991).  Resonant  cavity  applications  are  limited  to  small  volumes  of  material, 
which  could  present  a  problem  with  trying  to  measure  the  homogeneous 
response  of  substances  like  sands,  w-hich  are  highly  nonhomogcneous  at  small 
scales.  Waveguide  measurements  are  very  sensitive  to  how  the  sections  are 
mated  (and  how  repeatable  are  those  connections).  While  they  involve  ade¬ 
quate  volumes  of  material,  sample  preparation  is  not  likely  to  produce  a  urn- 
form  deasity  for  granular  materials,  as  the  substance  must  be  poured  into  one 
end  and  tamped  to  fill  the  volume.  Dielectric  seals  would  be  required  to 
provide  a  plane  interface  and  must  be  accounted  for  properly.  The  literature 
also  indicates  that  any  phase  shifts  because  of  nonideal  calibrations  are  nonlin¬ 
ear  functions  of  frequency  and  must  be  accounted  for  during  processing  of  the 
data  Waveguards  also  possess  a  cutoff  frequency  below  which  plane  waves 
cannot  propagate.  As  is  true  for  waveguides,  the  usual  coaxial  devices  (head¬ 
less  air  lines)  suffer  from  the  problems  of  uniform  sample  preparation  and,  in 
fact,  are  best  suited  to  measuring  solid  cylinders  of  material  that  must  be 
precisely  machined  to  provide  a  snug  fit  within  the  measurement  cavity. 
However,  coaxial  lines  allow  plane  waves  to  propagate  at  all  frequencies. 
Given  ail  of  these  considerations,  a  unique  coaxial  sample  holder  was  designed 
and  fabricated  for  use  in  these  experiments  to  provide  a  practical  tradeoff 
among  sample  volume,  repeatability  of  connections,  and  the  ability  to  prepare 
uniform  samples. 


Measurement  System 

A  schematic  of  the  measurement  system  usea  to  conduct  these  experiments 
is  shown  in  Figure  1  The  electronic  equipment  that  was  available  is  not  that 
which  would  normally  be  used  for  making  material  electrical  propeny  mea¬ 
surements,  but  it  was  successfully  adapted  for  this  use. 
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Figure  1.  Experimental  measurement  system 

The  Hewlett-Packard  85 IOC  Vector  Network  Analyzer  (Hewlett-Packard 
1991)  is  a  laboratory  device  designed  to  measure  the  magnitude  and  phase 
characteristics  of  electronic  networks  and/or  components.  It  consists  of  a 
detector  unit  that  measures  analog  signals  in  up  to  four  different  channels, 
down  converts  the  signals  to  an  intermediate  frequency  of  100  kHz,  performs 
an  analog-to-digital  conversion,  and  sends  those  data  to  a  data  processing  unit 
that  provides  a  display  of  the  desired  quantities.  The  HP  85 IOC  directs  the 
signal  source  (in  this  case,  an  HP  8340B  Synthesized  Sweeper)  to  generate  a 
sinusoidal  signal  at  a  specific  frequency  and  directs  a  te.st  set  (an  HP  8511 A 
Frequency  Converter)  to  act  as  the  receiver  for  up  to  four  channels  of  infor¬ 
mation  that  undergo  a  first  frequency  conversion  to  about  20  MHz  through  the 
process  of  harmonic  mixing  before  being  sent  to  the  detector  unit  described 
abo’/e. 

The  combined  dual-directional  coupler  and  frequency  converter  serves  as  a 
substitute  for  a  special  test  set  called  the  S-Parameter  Test  Set,  which  is 
usually  used  for  such  measurements  (Hewlett-Packard  1985).  The  S-Parame- 
ter  Test  Set  contains  power  splitters,  coupling  devices,  and  means  of  establish¬ 
ing  equal  circuit  electrical  lengths  where  appropriate  that  are  transparent  to  the 
user.  An  S-Paramcter  Test  Set  would  make  measurements  of  signals  sent 
through  the  transmission  ceil  in  both  directions.  As  long  as  one  believes  that 
the  material  response  is  not  dependent  on  which  direction  current  flows 
through  the  sample,  then  the  reverse  measurements  should  be  redundant. 
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I'igure  2  is  a  photograph  of  one  of  the  two  nearly  identical  sample  holders 
that  were  fabricated  for  these  measurements.  It  was  made  from  4.5-mm  brass 
stock  and  provides  a  stiuare  cross-sectional  coaxial  test  volume  that  is  7.5 
mm'  and  40  mm  long.  A  piece  of  copper  bus  wire  was  soldered  into  place  as 
the  center  conductor,  connecting  the  two  50-ohm  SMA-type  connectors  at 
each  end  of  the  sample  holder.  The  value  of  sample  volume  for  each  holder 
was  estimated,  by  length  measurements,  tt)  be  2.21  cm  . 


Figure  2.  One  of  the  brass,  coaxial  sample  holders 

One  area  of  concern  with  this  sample  holder  deals  with  the  development  of 
higher  order  transverse  modes.  Theory  says  that  a  coa.xial  device  should 
support  only  the  transverse  electromagnetic  mode  as  long  as  the  wavelength 
in  the  material  is  greater  than  twice  the  cavity  width  (Stratton  1941).  Thus, 
if  the  cavity  is  empty,  the  highest  frequency  allowed  for  these  measurements 
would  be 


3  -v  10’*’  t’w/scc 
v/T  .V  1.5  cmlcycle 


=  20  GHz 


(1) 


where 

V  =  wave  velocity  in  the  material 
X  =  wavelength 

«  =  inde.x  of  refraction  for  the  material  in  the  cavit> 
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c  ~  speed  of  light  in  a  vacuum 


while  for  pure  water 


3  X  10’° 
^  X  1.5 


=  2.24  GHz 


(2) 


Plane  wave  conditions  for  measurements  of  moist  soils  should,  therefore, 
exist  at  some  intermediate  frequency. 

A  Caron  Model  2065  Circulating  Bath  was  used  to  provide  sample  tem¬ 
perature  control.  The  bath  was  connected  to  a  hollowed  aluminum  heat 
source/sink  with  insulated  plastic  tubing.  As  the  circulating  50/50  antifreeze 
solution  controlled  the  physical  temperature  of  the  source/sink,  it  also  con¬ 
trolled  the  temperature  of  the  brass  sample  holder,  which  was  positioned  in  a 
slot  in  the  top  of  the  source/sink  and  covered  with  an  insulating  styrofoam 
cap. 


Governing  Equation 

Consider  the  geometry  of  the  coaxial  measurement  device  and  the  vol¬ 
tages  P;  in  all  regions  as  shown  schematically  in  Figure  3.  A  relatively  sim¬ 
ple  governing  equation  that  relates  the  sample  properties  and  dimension  to 
quantities  that  can  be  measured  by  the  network  analyzer  system  can  be 
developed  by  proper  consideration  of  the  boundary  conditions. 


^  X 

OUTER  CONDUCTOR 

1  I' 

IWT 

2*  2- 

3* 

CENTER  CONDUCTOR 

CABLE  .  Zp 

BB9BI 

''j  '  'l  CABLE  .  Zq 

H -  — H 

Figure  3.  Sample  geometry  and  voltage  and  current  notation 

Using  the  customary  notation  for  spatial  variation  of  transverse  voltages 
within  the  cables  and  sample  and  considering  a  snapshot  in  time  (without 
loss  of  generality,  let  time  be  zero),  one  can  write 
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X  ''  0 


V  =  Fi*  e  +  Kj 


-ill 


0  <  jc  <  J 


(3) 


;r  >  ^/ 


where  the  k's  are  the  complex  propagation  constants  and  contain  all  of  the 
material  electrical  property  information  (see  Appendix  A).  Similar  expres¬ 
sions  for  the  center  conductor  current  can  be  written  by  taking  advantage  of 
the  definition  of  the  characteristic  impedance  of  the  coaxial  line  as  being  the 
ratio  of  transverse  voltage  to  the  conduction  current  along  the  center  conduc¬ 
tor  (Stratton  1941). 


/j  =  J_  -  Kj'e 0  <x<d 


(4) 


A 


x>  d 


where  the  Z’s  are  the  characteristic  impedances  within  each  section  of  the 
line. 

The  boundary  conditions  that  can  be  applied  to  the  geometry  shown  in 
Figure  6  are  that  the  tangential  field  components  at  the  interfaces  are  contin¬ 
uous  and,  from  the  steady-state  continuity  relationship,  that  the  currents  are 
also  continuous  across  the  interfaces.  In  equation  form,  these  boundary  con¬ 
ditions  reduce  to 
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V^,  =  @  ^  =  0  ; 


(5) 


=  V,@x  =  d  : 


(6) 


/,  =  /j  @  x  =  0  : 


-  v','  _  -V2  -  v; 


2o 


(7) 


I2  =  ®  X  =  d  : 


-  v-“>^  .v; 


(8) 


Not  being  able  to  measure  V*'2  or  V'2  but  certainly  being  able  to  obtain  a 
reading  on  the  other  voltages,  the  two  voltage  amplitudes  inside  the  material 
can  be  eliminated  by  adding  together  Equations  S  and  7  and  by  taking  their 
difference; 


v:  =  i 


'W  "  ''i'  *  :|('W  -  '^r) 


(9) 


v;  -  4. 


%  -  yi  -  ^{Virp  -  V,') 


(10) 


Substituting  these  results  into  Equation  6,  one  has 


v;  =  1 
^  2 


*ik/l 


(11) 


1 

2 


Zs 


Vinp  "  I'W  -  ''l, 


(iw  -  ni 


-ik/I 
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Normalizing  lo  V,„|,  and  collecting  terms. 


i-t  — 

y 

tnp 

cos  ^k^cfj  +  / 

~ 

tnp 

V 

mp 


(12) 


A  similar  relationship  can  be  written  from  Equation  8 


cos  (k^d)  +  /• 


sin 


2o  Kn^ 


(13) 


It  is  further  convenient  to  eliminate  the  characteristic  impedances  by  solvi;,g 
each  of  the  last  two  equations  for  and  setting  the  two  expressions  equal: 


v; 

' 

cos  {k^d^  +1 

'*-1 

V 

sin(iff) 

r  1 

V 

imp 

sin  {k,<i)  - 

^  imp 

v; 

~  T~ 

ifip 

cos  (k^dj 

(14) 


Simplifying, 


cos  (*  rf) 


(15) 


Now,  following  an  appropriate  calibration  of  the  measurement  apparatus, 
the  network  analyzer  effectively  measures  V,A^(,p  and  They  are 

designated,  using  standard  ”8",  or  scattering  parameter  notation,  as  S,,  and 
Sji,  respectively.  The  definition  of  these  "S"  parameters  requires  that  there 
be  no  current  traveling  into  the  sample  from  the  right  side  of  Figure  3.  This 
implies  that  all  circuitry  and  connectors  on  the  transmission  side  of  the  sam¬ 
ple  are  perfectly  matched  so  that  no  energy  is  reflected  back  to  the  sample. 
Although  this  is  an  assumption  for  this  derivation,  the  fact  that  proper  equip¬ 
ment  was  not  available  for  these  measurements  and  that  a  less  than  desirable 
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calibration  technique  was  used  (described  in  the  next  section)  makes  this 
idealized  configuration  quite  unlikely. 

Using  the  relationships  between  the  propagation  constant  and  material 
properties  (Appendix  A)  and  assuming  nonmagnetic  materials,  one  can  then 
write  dowTi  the  final  form  of  the  governing  equation: 


cos 


2S, 


=  fi' 


(16) 


This  form  of  the  governing  equation  was  reported  by  Campbell  (1990).  An 
earlier  derivation,  although  arriving  at  a  different  solution  technique,  has  also 
been  recently  discovered  (Hewlett-Packard  1985). 

The  process  for  calculating  values  of  the  complex  dielectric  constant  can 
now  be  outlined; 

a.  At  a  given  frequency,  measure  the  complex  values  of  S,,  and  Sj,. 

b.  Compute  cos"'  (H’)  =  aid  —  . 

'  c 

c.  Square  both  sides  of  this  expression  and  solve  for  e. 

The  procedure  is  not  as  straightforward  as  it  seems  because  of  a  problem 
with  principal  values  of  the  cos  '  function.  One  method  ‘or  solving  for  the 
cos  '  function  yields  the  following  (Zucker  1972): 


cos"'  =  -i  In  (re'^)  =  6  -  i  In  (r) 


where 


rc*  =  -  1 

A  great  deal  of  time  was  spent  developing  a  calculation  procedure  for  Equa¬ 
tion  17  that  would  yielJ  the  correct  values  of  e  and  e*.  However,  it  was 
later  discovered  that  the  calculation  .outi  le  for  the  inverse  complex  cosine 
function  that  resides  in  the  HP  BASIC  compiler  within  the  HP9000  system 
controller  gave  exactly  the  same  results.  The  only  difficulty  is  that,  depend¬ 
ing  on  where  lies  in  complex  space,  the  cos  '(W^)  value  computed  using 
the  available  routines  •etums  real  and  imaginary  components  that  are  both 
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negative,  or  having  opposite  signs.  A  set  of  rules  had  to  be  developed  to 
select  the  sensible  results;  namely,  those  that  would  yield  positive  values  of  e 
and  €  . 

The  real  part  of  udy^/c  is  a  measure  of  the  number  of  wavelengths 
(expressed  in  radians)  of  signal  that  can  be  contained  w-ithin  the  sample  holder 
at  the  given  radial  frequency.  The  value  of  6  that  is  returned  by  the  above 
solution  algorithm  is  a  principal  value,  being  bounded  by  ±ir.  However,  any 
multiple  of  can  be  added  to  6  and  give  the  same  result.  It  rums  out  that 
there  is  a  sensible  way  of  determining  how  many  multiples  of  2n  to  add.  One 
can  either  estimate  the  properties  of  the  material  to  establish  the  correct  value 
of  6,  or  one  can  use  the  value  of  "group  delay,"  which  is  calculated  by  the 
network  analyzer  to  establish  the  correct  value  of  6.  Group  delay  is  defined 
as  the  signal  transit  time  through  a  test  device  (Hewlett-Packard  1991).  WTiat 
the  designers  of  the  network  analyzer  assumed  is  that  the  phase  for  the  wave 
in  the  material  is  a  term  like  "wt".  Therefore,  group  delay,  or  the  derivative 
of  phase  with  respect  to  radial  frequency  is  a  number  with  dimensions  of  time 
and  interpreted  as  the  time  required  to  transit  the  sample.  An  estimate  of  the 
initial  value  of  the  number  of  cycles  (multiples  of  27r)  to  add  to  6  was  taken  as 
the  truncated  product  of  group  delay  and  frequency. 

Appendix  B  contains  a  listing  of  the  HP  BASIC  code  tliat  was  developed 
^v  the  author  to  control  the  data  collection  system,  create  a  permanent  data 
file  of  amplitude  and  phase  measurements  as  a  function  of  frequency,  and 
produce  a  real-time  plot  of  computed  complex  dielectric  constant  values, 

Experimental  Procedure 

Calibration 

Because  of  the  arrangement  and  type  of  measurement  apparatus  used  for 
these  studies,  the  most  appropriate  calibration  was  found  to  be  the  One-Path, 
2-Port  calibration  combined  with  the  adapter-switching  technique  (Hewlett- 
Packard  1991).  During  this  procedure,  reflection  and  transmission  cables 
were  first  connected  with  a  thru  adapter  from  an  HP85052B,  3.5inm 
Calibration  Kit.  and  forward  transmission  and  forward  match  measurements 
were  conducted.  Next,  an  equivalent  adapter  (male-female)  from  the  same  Kit 
was  connected  to  the  input  (reflection)  cable,  and  three  reflection  standards 
were  measured.  These  standards  were  an  offset  short,  an  offset  open,  and  a 
broadband  load,  all  from  the  same  Calibration  Kit.  The  net  result  of  the 
calibration  procedure  is  that  the  Sji  measurement  effectively  takes  place  at  the 
bead  interface  in  the  SMA  connector  on  the  reflection  side  of  the  sample 
holder,  and  the  S21  measurement  effectively  takes  place  at  the  bead  interface 
in  the  connector  on  the  transnussion  side  of  the  holder. 

Because  of  the  physical  separation  of  the  calibration  planes  and  the  inter¬ 
face  between  the  sample  and  the  connector  on  the  end  of  the  holder,  an 
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Because  of  the  physical  separation  of  the  calibration  planes  and  the  inter¬ 
face  between  the  sample  and  the  connector  on  the  end  of  the  lioluer,  an 
additional  phase  shift  needed  to  be  added  to  the  c?,,  parameter  and  subtracted 
from  the  S,^  parameter  to  make  the  measured  quantities  consistent  with  what 
was  defined  in  the  earlier  derivation  of  the  governing  equation.  Referring  to 
the  sketch  in  Figure  4,  the  reflected  signal  must  travel  an  additional  distance 
of  2a  beyond  the  calibration  plane  to  be  recorded- in  the  frequency  converter. 
The  transmitted  signal  travels  an  additional  distance  of  a-<-b.  I  or  these  mea¬ 
surements,  a  =  b,  and  the  required  phase  shift  was  found  to  be  about  80  pico¬ 
seconds  from  measurements  of  the  empty  holder. 


Figure  4.  Calibration  and  measurement  planes 


Sample  preparation  and  measurements 

Test  samples  can  be  prepared  from  any  reasonably  fine-grained  soil  or 
from  any  liquid  that  does  not  chemically  interact  with  the  sample  holder 
and/or  center  conductor.  Each  soil  sample  should  normally  be  made  as 
homogeneous  as  possible  by  drying,  pulverizing,  and  mixing  thoroughly. 
Samples  should  then  be  taken  to  determine  the  usual  engineering  properties 
of  grain-size  distribution  and  plasticity  indices. 

The  normal  measurement  procedure  for  solid/liquid  mixtures  is  the 
following: 

a.  Lower  the  temperature  of  the  heat  source/sink  to  the  lowest  value 
desired. 

b.  Carefully  weigh  the  dry,  empty  sample  holder. 
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c.  The  sample  holder  should  be  filled  with  a  nearly  saturated  sample  of 
the  mixture  that  has  been  prepared  m  an  appropriate  container  (for 
soil  water  mixtures,  resealable  plastic  bags  were  found  to  be  conven¬ 
ient).  The  holder  can  be  tapped  against  a  flat,  hard  surface  to  uni- 
foiTnly  settle  the  sample. 

d.  The  sample  and  holder  should  be  carefulh  weighed,  then  sealed  and 
connected  to  the  measurement  system. 

e.  The  measurement  program  containing  the  dielectric  property  calcula¬ 
tion  algorithm  should  then  be  executed,  which  will  cause  the 
frequency  range  for  the  dual  directional  coupler  being  used  to  be 
stepped  through  and  calculations  of  e  and  e  made. 

/.  The  temperature  of  the  sample  can  then  be  increased  and  the  mea¬ 
surement  program  executed  again.  For  each  new  temperature  setting, 
the  sample  should  be  given  at  least  10  min  to  reach  equilibrium  once 
the  bath  has  reached  the  desired  setting.  The  10-min  wait  was  found 
experimentally  to  be  more  than  adequate  for  sample  equilibrium  to  be 
achieved. 

g.  At  the  end  of  any  particular  measurement,  the  holder  cover  should  be 
loosened  to  provide  a  better  opportunity  for  the  sample  to  dry  out. 

The  sample  and  holder  should  be  periodically  weighed  to  detennine 
approximate  values  of  moisture  content.  When  the  sample  is  thought 
to  be  at  about  the  desired  volumetric  moi.sture  content,  the  holder 
should  be  resealed  and  reconnected,  the  sample  temperature  should  be 
taken  back  to  the  lowest  desired  temperature,  and  the  whole  measure¬ 
ment  sequence  should  be  repeated. 

h.  .After  all  measurements  have  been  conducted,  the  sample  should  be 
carefully  removed,  dried,  and  weighed  to  determine  its  original  dry 
density  and  to  facilitate  calculations  of  volumetric  moisture. 
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3  Sanity  Checks 


To  develop  confidence  in  the  experimental  procedure,  materials  were  mea¬ 
sured  whose  responses  are  known  experimentally  or  can  be  modeled  very 
accurately.  These  measurements  also  served  to  confirm  the  amount  of  phase 
shift  required  to  collapse  the  calibration  and  measurement  planes  as  discussed 
in  the  previous  chapter.  Data  were  collected  for  empty  sample  holders,  sam¬ 
ples  of  distilled  water  and  ethylene  glycol  (ethanediol),  and  a  sample  of  silty 
soil.  Each  is  discussed  in  the  following  sections. 


Empty  Holder 

The  first  question  to  answer  is  'How  well  docs  this  cxpciimcntal  procedure 
measure  the  properties  of  air?"  Ideally,  the  real  pan  of  the  complex  dielectric 
constant  for  dry  air  should  be  unity,  and  the  imaginary  pan  should  be  zero. 
Figure  5  shows  the  results  of  measuring  empty  sample  holders.  Considering 
all  of  the  experimental  variables  that  could  have  an  impact  on  any  measure¬ 
ment  such  as  measurement  system  drift,  connector  and  connection  variability, 
outer  conductor  discontinuities  created  by  the  sample  holder  cover,  etc.,  the 
agreement  of  entity  bolder  (air)  measurements  with  theoretical  predictions  is 
quite  remarkable. 

The  venical  dashed  line  at  2  GHz  separates  measurements  made  in  two  dif¬ 
ferent  frequency  ranges.  A  single  dual  directional  coupler  (Figure  1)  that 
covers  the  entire  50  MHz  -  26.5  GHz  frequency  range  that  is  available  from 
the  source  does  not  exist.  Two  different  couplers  were  used  whose  ranges  of 
operation  overlappied  at  2  GHz. 

The  anomalous  behavior  centered  about  3.75  GHz  and  7.5  GHz  is  asso¬ 
ciated  with  the  mathematical  solution  technique.  At  3.75  GHz,  the  empty 
sample  holder  should  be  spanned  by  about  a  half  of  a  wavelength;  at  7.5  GHz, 
it  should  hold  about  a  full  wavelength,  and  so  on.  The  number  of  wave¬ 
lengths  contained  within  the  sample  holder  for  a  material  with  no  loss  is  pro¬ 
portional  to  (<  dv^cor  the  cosine  function  argument  in  Equation  16.  A  half 
wavelength  in  the  sample  holder  will  result  in  the  complex  vector  representing 
the  right-hand  side  of  Equation  16  approaching  the  negative  real  axis  in  com¬ 
plex  space.  A  full  wavelength  corresponds  to  the  vector  having  made  a 
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complete  revolution  in  complex  space  and  approaching  the  positive  real  axis 
But  these  two  conditions  are  branch  cuts  for  the  cos'^(z}  function  and  are. 
therefore,  inherently  uastable.  It  has  been  recently  discovered  that  other 
researchers  have  experienced  this  anomalous  behavior  and  that  at  least  one 
alternative  solution  technique  has  been  successfully  developed  (Baker-Jarvis, 
Vanzura,  and  Kissick  1990).  Future  improvements  to  the  measuiement  and 
calculation  scheme  reponed  herein  will  probably  include  one  or  more  of  these 
alternative  solutions. 

Much  of  the  data  that  follow  in  this  chapter  show  the  same  anomalous 
behavior  whenever  'rteger  multiples  of  half  wavelengths  exist  within  the  sam¬ 
ple  holders.  The  severity-  of  these  peaks  decreases  with  increasing  frequency, 
because  at  low  frequencies  the  computation  of  t  from  Equation  16  involves  a 
division  by  very  small  numbers.  Any  inaccuracies  are  amplified  at  low  fre¬ 
quencies  (<  1  GHz).  This  fact  also  had  a  negative  impact  on  low-frequency- 
measurement  calibrations.  Instrumenution  drift  was  often  enough  to  ruin  a 
calibration  at  the  lowest  frequencies  in  just  a  few  hours  Better  measurement 
equipment  (namely,  an  S-parameter  test  set  and  a  high-quality  set  of  cables 
and  connectors)  would  help  alleviate  this  low-frequency  measurement  pro¬ 
blem 


Water 

Measured  values  of  the  complex  dielectric  constant  for  samples  of  distilled, 
deionized  water  are  shown  in  Figure  6.  For  comparison,  Ray’s  (1972)  empir¬ 
ical  fits  to  a  large  amount  of  pure  and  freshwater  dau  are  duplicated  here. 
Several  observations  are  in  order.  First  of  all,  the  comparison  of  measured 
data  to  the  accepted  standard  is  quite  good  (less  then  4  percent  difference), 
especially-  for  the  imaginar,.  term.  The  sample  holder  volume  was  evidently 
not  completely  filled  with  water  because  of  the  fact  that  the  holder  cover  must 
be  pressed  down  upon  the  water  meniscus,  causing  water  to  be  squeezed  from 
the  holder  in  an  unpredicuble  way.  It  is  also  possible  that  some  evaporation 
of  the  sample  took  place  during  the  2-  to  3-hr  measurement  period.  Both  of 
these  phenomena  could  account  for  t’ne  lower  measured  values  of  c  .  The 
hypothesis  that  the  sample  holder  was  never  completely  full  was  supponed  by¬ 
weight  measurements  that  indicate  the  holders  were  never  more  than  95  per¬ 
cent  full  during  these  tests. 

Another  observation  is  that  because  of  the  good  correlation  with  idealized 
data,  the  impact  of  higher  order  modes  must  be  relatively  small  as  one  w-ould 
expect  (Stratton  1941).  As  pointed  out  at  Equation  2,  higher  order  transverse 
modes  should  appear  above  2.3  GHz  for  pure  water,  but  these  data  thaf  extend 
to  4  to  5  GHz  continue  to  behave  well,  f'^'ater  data  are  cut  off  at  4  to  5  GHz 
because  the  noise  floor  of  the  network  analyzer  measurement  system  w-as 
reached.)  Tlie  same  anomalies  observed  at  half-wavelength  intervals  in  the 
emptv  holders  also  occur  in  the  w-ater  with  very  little  impact  being  observed 
beyond  one  and  a  half  wavelengths. 
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Figure  5.  Empty  sample  holder  measurements 

The  low-frequency  tail  on  e'  may  possibly  be  associated  with  the  ionic  con¬ 
ductivity  loss  mechanism  as  a  result  of  either  trace  quantities  c  f  ions  that 
escape  the  deionization  process  cr  as  a  result  of  ions  that  originate  from 
unclean  and/or  oxydized  outer  or  inner  conductor  surfaces.  The  magnitude  of 
low-frequency  losses  are  approximately  those  associated  with  freshwater  con¬ 
ductivity,  but  the  actual  loss  mechanism  has  to  be  viewed  as  purely 
speculative. 


Ethylene  Glycol 

Data  were  collected  on  several  different  alcohols:  ethylene  glycol 
(ethanediol),  methanol,  and  ethanol.  However,  because  of  surface  tension 
problems  md  rapid  evaporation  rates,  only  the  results  for  ethylene  glycol 
(HOCHjCHjOh,  molecular  weight  62.07,  specific  gravity  1.1088  at  20  “C 
relative  to  water  at  4  “C,  <  .05  percent  water)  are  reproduced  here  (Figure  7). 
Measured  values  of  both  f  and  f"  are  consistently  lower  than  .Iordan’s 
(Jordan,  Sheppard,  and  Szwamowski  1978)  results  for  frequencies  greater 
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Hgure  6  Dispersion  curves  for  water 
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Figure  7.  Dispersion  measurements  of  ethylene  glycol  (Continued) 
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than  1  GHz.  This  is  not  a  surprising  observation  in  that  the  experimental 
procedure  does  suffer  from  the  shortcoming  that  it  is  difficult  to  completely 
fill  the  sample  holder  volume  with  liquids  because  of  surface  tension  prob¬ 
lems,  and  it  is  difficult  to  maintain  the  volume  of  liquid  during  the  2-  to  3-hr 
measurement  period  because  of  evaporation.  Although  ethanediol  was  least 
affected  by  these  problems,  sample  weights  for  the  high-frequency  measure¬ 
ments  indicate  that  the  holder  was  90  to  95  percent  full  during  the  experi¬ 
ments.  At  low  frequencies,  measurements  revealed  another  loss  mechanism, 
most  likely  ionic  conductivity. 

In  summary,  results  of  measurements  of  air,  water,  and  ethylene  glycol  are 
quite  good,  resulting  in  a  high  degree  of  confidence  in  the  experimental  proce¬ 
dure.  Differences  between  these  measurements  and  more  precise  measure 
ments  reported  in  the  literature  can  be  explained  by  physical  and'ot  chemical 
factors.  One  must  conclude,  however,  that  improvements  would  have  to  be 
made  to  the  experimental  procedure  before  it  could  be  used  to  accurately 
measure  the  electrical  properties  of  a  variety  of  fluids.  In  particular,  it  would 
be  necessary  to  find  a  way  to  completely  fill  the  sample  holder  and  to  main¬ 
tain  that  fluid  level  throughout  the  measurements.  The  only  current  concerns 
with  measurements  oi  heterogeneous  mixtures  of  soil  and  water  are  cleanliness 
of  the  sample  holder  and  uneven  dr>  ing  of  the  sample  throughout  its  volume. 


A  Silty  Soil 


As  an  example  of  how  data  can  look  for  heterogeneous  mixtures  such  as  a 
moist  soil,  consider  the  data  shown  on  Figure  8.  These  data,  representing  two 
volumetric  moisture  content  extremes  (11.7  and  43.8  percent),  are  from  mea¬ 
surements  made  at  20  "C  on  a  sample  of  very  clean  and  somewhat  poorly 
graded  silt.  More  information  on  the  properties  of  this  soil  will  be  available 
in  the  third  report  of  this  series. 

The  results  of  these  measuremems  are  certainly  sensible  For  example,  the 
real  pan  of  the  complex  dielectric  constant  should  be  less  than  that  of  pure 
water,  as  shown  in  the  previous  figure.  At  the  same  time,  the  loss  term 
should  contain  evidence  of  the  relaxation  peak  beyond  10  GHz,  which  it  does, 
and  should  reveal  ?  much  greater  direct  current  conductivity  contribution 
below  1  GHz  because  of  the  soil  chemistry  and  the  fornration  of  free  ions  in 
the  soil /water  mixture. 
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4  Summary 


The  response  of  moist  soils  (or  any  other  moist  heterogeneous  mixture)  to 
active  microwave  sensors  is  controlled  by  surface  geometry  and  the  dielectric 
properties  of  the  medium.  A  review  of  the  literature  reveals  that  what  data  do 
exist  on  dielectric  properties  suffer  from  a  number  of  shoncomings.  In  most 
cases,  the  soils  are  simply  not  properly  characterized  in  terms  of  their  physical 
parameters  such  as  the  distribution  of  particle  sizes,  the  dry  density  under 
which  tests  were  conducted,  or  the  temperature  of  the  soil,  all  of  which  have 
some  impact  on  the  dielectric  properties.  The  volumetric  moisture  content  of 
the  soil,  which  is  usually  the  most  important  factor  controlling  dielectric 
response,  is  often  not  specified  or  is  not  computable.  And.  finally,  the  data 
are  usually  collected  at  a  limited  numbe*’  of  frequency  values. 

In  an  effort  to  supplement  those  data  on  the  complex  dielectric  properties 
of  moist  soils  that  have  been  adequately  documented,  a  coaxial  measurement 
apparatus  was  fabricated  to  allow  electrical  property  measurements  in  nom¬ 
inally  moist  soils  to  frequencies  above  10  GHz.  Sample  temperatures  were 
controlled  by  an  external  bath.  The  algorithm  used  to  calculate  the  dispersive 
complex  dielectric  constant  was  verified  by  measuring  the  response  of  liquids 
whose  properties  are  well  documented  and  the  response  of  air  (empty  sample 
holder). 
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Appendix  A 

Skin  Depth  Calculations 


Obviously,  skin  depth  effects  preclude  the  utility  of  practical  elec¬ 
tromagnetic  devices  to  penetrate  to  distances  comparable  vvitli  acoustic  devices, 
but  nevertheless,  low-frequency  devices  have  some  utility  (Berlin  ct  al.  1986; 
Farr  et  al.  1986;  Schaber  ct  al.  1986).'  As  a  means  of  quantifying  the  pene¬ 
tration  of  electromagnetic  waves  into  earth  media  in  terms  of  the  complex 
dielectric  constant  of  the  media,  one  should  consider  the  following  arguments. 
The  onc-dimcnsional  propagation  of  an  electromagnetic  wave  into  some 
medium  can  be  described  by  an  amplitude  function  such  as 


e‘ 


(it*  -v*) 


(Al) 


where 

i  =  symbol  designating  the  imaginary  quanuty 
k  =  complex  wave  number 
X  =  distance 
0)  =  radial  frequency 
t  =  lime 
Furthermore,  let 

*  i  ^  (A2) 

2 


'  Refem.ces  cited  in  this  appendix  are  located  a:  the  end  of  the  main  text. 
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where 

P  =  pha<ic  consiani  and  governs  how  a  lossless  medium  propagates  a 
wave 

a  =  aitenuauon  factor  and  controls  how  the  amplitude  of  the  wave  de¬ 
creases  wiLh  distance  traveled  through  the  medium 


TTicn  the  amplitude  of  the  wave  is 


(A3) 


i.c.,  a  traveling  wave  combined  with  an  attenuation  factor. 

If  skin  depth  is  defined  as  that  distance  over  which  the  wave  amplitude 
decreases  by  a  factor  of  l/c,  then  the  skin  depth  8  is 

8  H  i-  (A4) 

a 


But  the  wave  number  can  also  be  wntten 


0)  _  toA' 
V  c 


where 

V  =  phase  velocity  in  the  medium 

c  =  speed  of  light 

N  =  complex  index  of  refraction 
=  n  -t-  in’ 

Thus 


(A5) 


8  =  2=  (A6) 

a  cun" 

or  the  skin  depth  is  inversely  proportional  to  the  imaginary  part  of  the  complex 
index  of  refraction. 

However,  this  study  focuses  on  the  complex  dielectric  constant  e  as  the 
measure  of  soil  electrical  charactenstics 


A2 
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e  =  e'  *  iz"  =  A'‘ 

Following  several  substitutions  and  simplifications,  one  can  show  that 


(A8) 


where  /vfl  is  the  free-spacc  wavelength.  A  visualization  of  this  equation  for  the 
range  of  values  of  permittivity  z  and  loss  tangent  z  Iz  normally  found  in  soils 
is  provided  in  Figure  Al.  This  figure  does  not  represent  any  new  science,  as 
the  same  information  has  been  prc.sentcd  in  other  forms  (von  Hippcl  1954).  It 
is,  however,  a  more  condensed  way  of  visualizing  the  skin  depth  that  has  not 
been  observ'cd  elsewhere  in  the  literature.  Figure  Al  is  valid  for  any  interpre¬ 
tation  of  loss  tangent;  i.e.,  whether  one  is  speaking  of  a  loss  mechanism 
because  of  actual  charge  migration  or  one  because  of  dielecLhc  relaxation. 

Another  useful  skin  depth  relationship  found  m  the  literature  is  the  nomo¬ 
graph  generated  by  Albrecht  0966)  and  reproduced  on  Figure  A2.  Albrecht 
first  developed  an  empirical  relationship  that  associated  ground  conductivity  to 
gravimetric  moisture  content  and  ground  temperature  at  low  frequencies.  He 
then  related  skm  depth  to  wavelength  and  conductivity  by  the  usual  equation 
for  highly  conductive  materials  (assumes  that  the  conductivity  phenomenon 
dominates  the  loss  mechanism  in  soils  at  low  frequencies); 


(A9) 

npco 

where  the  units  used  are  MKS,  and  the  skin  depth  is  in  meters  and  tlie  ground 
conductivity  is  in  units  of  mhoAneter,  The  magnetic  permeability  was  taken  to 
be  47cX1{7^  Henrys./meter.  The  way  to  use  this  nomograph  is  to  select  a  mois¬ 
ture  content  and  move  on  a  horizontal  line  to  the  appropriate  temperature 
curve.  Then  m.ove  on  a  vertical  line  to  the  appropriate  frequency  curve.  From 
that  point,  a  horizontal  line  to  the  right-hand  scale  reveals  the  desired  value  of 
skin  depth.  There  is  no  way  to  relate  Albrecht’s  results  to  those  of  Figure  Al, 
because  it  begins  with  an  unknown  empirical  relationship  among  moisture, 
temperature,  and  conductivity. 


Hoekstra  and  Delaney  (1974)  also  pubbshed  some  useful  data  on  plane 
wave  attenuation  in  moist  soils  at  higher  frequencies.  These  data  are  shown  in 
Figure  A3.  'Iliey  concluded  that  attenuation  was  relatively  independent  of  soil 
fabric  and  that  either  passive  or  active  microwave  sensors  could  obtain  a  mea¬ 
surable  response  from  only  about  the  first  .5  cm  of  ground  below  the  surface. 
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Figure  A2  Skin  depth  nomograph  (from  Albrecht  (1966 
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Appendix  B 

Program  Listing  and  Sample 
Output 


The  following  program,  written  in  the  HP  BASIC  language  and  executed  in 
an  interactive  mode  on  an  HP9000  Series  360  system  controller,  drives  the 
vector  network  analyzer-based  system  used  to  collect  complex  reflection/ 
traitsmission  parameters,  store  those  data  on  permanent  files,  and  generate  a 
real-time  hard  copy  output  of  calculated  complex  dielectric  constant  values. 

An  example  output  sheet  is  also  included. 
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!  PROGRAM  NAME;  EPSILONS 


10 
11 

12  !  A  multipurpose  CODE  FOR  1)  COLLECTING  COAXIAL  REFLECTION  AND 

13  !  TRANSMISSION  DATA  AND  CALCULATING  BOTH  THE  REAL  AND  IMAGINARY 

14  !  PARTS  OF  THE  COMPLEX  DIELECTRIC  CONSTANT,  OR  2)  DISPLAYING 

15  1  THE  RESULTS  OF  PREVIOUSLY  COLLECTED  DATA 

16  I 

17  1  THE  HIGH  FREQUENCY  OPTION  USES  AN  HP11692D  DIRECTIONAL  COUPLER 

18  J  AND  COVERS  THE  1.5-17.5  GHz  RANGE. 

19  I 

20  !  THE  LOW  FREQUENCY  OPTION  USES  AN  HP778D  DIRECTIONAL  COUPLER 

21  1  AND  COVERS  THE  45  MHz  TO  2.045  GHz  RANGE. 

23  ! 

24  1  USES  THE  HP9000  ACS  ALGORITHM  WITH  SIGN  ADJUSTMENTS 

25  I 

26  OPTION  BASE  1 

30  COMPLEX  Sll( 1 : 801)  ,S21( 1:801) ,Rh8,LhB,Rt,Eps 

40  INTEGER  Preaicble , Bytes , Phflcorr 

41  REAL  Blnary( 1:801, 1 : 2) 

50  DIM  TtlSieO) ,SeclineS(60) ,EpBp( 1:601) ,Ep8pp( 1:801) 

60  ASSIGN  @Nwa  TO  716 

70  ASSIGN  eNw«_data  TO  7 16, -FORMAT  OFF 

00  LOCAL  @Nwa 

90  PRINTER  IS  CRT 

100  GINIT 

110  PLOTTER  IS  CRT, "INTERNAL" 

12C  GRAPHICS  ON 

130  FdS=DATES(TIM£DATE; 

131  FtS'TIMES(TIM£DATE) 

150  CLEAR  SCREEN 

160  PRINT  FdS,Ft$  )  D.ATE  AND  TI.ME  AT  HECINNINC 

OF  MEA£UR£,MENT 

170  BEEP  4000, .2 

171  D*=4 

173  PRINT 

174  PRINT  "INSTRUMENT  STATE  1  IS  ASSOCIATED  WITH  THE  HIGH  FREQUENCY  COUPLER" 

175  PRINT  "INSTRUMENT  STATE  2  IS  ASSOCIATED  WITH  THE  LOW  FREQUENCY  COUPLER" 

176  PRINT 

177  INPUT  "WHICH  INSTRUMENT  STATE  DO  YOU  WANT  TO  USE?",Nl 

179  Pi«=4,5E*7 

180  Ff=2.045Et9 

181  Np-401 

162  IF  Nl=2  THEN  GOTO  190 

183  F1-1.5E+9 

184  Ff«1.75E*10 

185  Np-801 

189  I 

190  ASSIGN  eHplb  TO  7 

191  REMOTE  eapib 

192  ABORT  @Hpib 

193  CLEAR  @Nwa 

194  CLEAR  717 

196  OUTPUT  @Nwa; "WAIT;RECA"iVALS(Nl)t";"  1  RECALL  INSTRUMENT  STAX 

E 

197  OUTPUT  @Nwa;  "WAIT;CHANl;WAIT;S12;LOCKAl;tIUMEB2,-DENOAl;KEDD; "  1  REDEFINE  S 

12  t  S22  BECAUSEA  1-PATH,  2-PORT  CAJ,  IS  BEING  USED 

198  OUTPUT  @Nwf.;  "WAIT;622;LOCKA1;HUMED1;DENOA1;R£DD;" 

199  OUTPUT  eNwii;  "WAIT;Sll;  " 

200  OUTPUT  eNwa;"CALS"tVAL5(Nl)i";"  I  RECALL  CAL  SET 
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2C1  OUTPUT  @Nwa;  "WAIT;  PRIKTYPECO;.R; "  !  PRINT  COLORS 

202  BEEP  4000, . 2 
20J  PRINT 

204  INPUT  "ARE  YOU  COLLECTING  NEW  DATA  (1)  OR  RESTORING  AN  EXISTING  FILE  (2)? 
,Ndf 

20b  IF  Ndt=l  THEN  GOTO  246 
206  : 

707  !  RETRIEVE  AN  EXISTING  DATA  FILE  FOR  FURTHER  ANALYSIS 

208  : 

209  BEEP  4000,. 2 

210  PRINT 

211  PRINT  "PREPARE  FOR  A  DUMMY  MEASUREMENT  TO  PROPERLY  RESET  THE  INSTRUMENT 
TATE" 

212  PRINT  "WHEN  YOU  ARE  READY,  PRESS  THE  'CONTINUE'  KEY" 

213  PAUSE 

214  Wtinie“60 

215  IF  Ni-=2  THEN  Wtlme  =  62 

216  OUTPUT  @Nwa; "SING; " 

217  WAIT  3 

218  TRIGGER  @Nwa 

219  WAIT  Wtime  ;  TIME  REQUIRED  FOP.  A  SWEEP  WITH  AVERAGING 

220  'TRIGGER  @Nwa 

221  WAIT  Wtlme 

222  BEEP  4000,-2 

223  PRINT 

224  INPUT  "WHAT  IS  TTIE  NAME  OF  THE  EXISTING  DATA  F ILE7 ”  , Nainef  S 

225  ASSIGN  en  TO  "D I SSDATA/ "iName'S ; FORMAT  OFF 

226  OUTPUT  eNwa;"CHANl;FORM3;lNFURAWl,-" 

227  ENTER  eFl;Preatnble;BytGB 

228  M«Byte8/16 

229  PRINT  "N'JMBER  OP  DATA  POINTS  tN  FILE=",K 

230  REDIM  Binary ( 1 :M , 1 : 2 ; 

231  REDIM  Sll( 1:M),S21(1:M) 

232  ENTER  @Fl;Binary( * ) 

233  OUTPUT  @Nwa_data;PreatRble,Eyte8,Binary( * ) 

238  OUTPUT  @thrfa;"WAlT;CHAN2;FOKMb;INPURAW2,'" 

239  ENTER  eFl;Preanible;Byte8;Blnary(  *  ) 

240  OUTPUT  gNwa  data; ProamblB; Bytep; Binary (* ) 

241  ASSIGN  eKl  TO  • 

242  GOTO  302 

243  1 

244  I  PREPARE  TO  MAKE  A  NEW  SET  OF  MEASUREMENTS 

245  1 

246  OUTPUT  @Nwa; "OPEP; " 

263  Naraef$-Fd$( l,2]iFdS(4,6)SFd$( 11, 11 )6FtS( 1 , 2 ) SFtS [ 4 , 5 ) 

264  PRINT 

265  PRINT  "THE  BINARY  DAT  FILE  NAME  IS  ",NajDGlS 

266  IF  Nl«l  THEN  GOTO  270 

267  -CREATE  BDAT  "/DIS3DATA/ "tNamof $ , 1 , 6440  J  CREATE  DATA  FILE 

268  GOTO  271 

270  CREATE  BOAT  "/DISSDATA/ "tNanaf S , 1 , 1 2848  I  CREATE  DATA  FILE 

271  ASSIGN  ®F1  TO  " /DI3SDATA/ "ANamotS; FORMAT  OFF 
2/3  WtllD«»80 

274  IF  Ni-2  THEN  Wtlire-62 

275  OUTPUT  «Nwa;"SING; " 

276  WAIT  3 

277  TRIGGER  @Nwa 

278  WAIT  Wtlme  1  TIME  REQUIRED  FOR  A  SWEEP  WITH  AVERAGING 

279  TRIGGER  @Nwa 

280  WAIT  Wtimo 
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281 

262 

1  STORE  A  NEW  DATA  FILE 

203 

1 

264 

OUTPUT  @Nwa;  "CHANl;  F0RM3  ;OUTPRAWl ; 

265 

ENTER  @Nwa_data;Preanible;Bytes 

28C 

M'Bytes/ 16 

287 

REDIM  Binary(l:M,l:2) 

266 

BEDIM  Sll(l'M) ,S21(1:M) 

289 

ENTER  @Nwa_data;Binary( * ) 

294 

OUTPUT  @Fl;Preaiiible;ByteB;Binary(  •  ) 

295 

OUTPUT  eNwa;"WAIT;CHAN2;FORM3;OUTPRAW2; " 

296 

ENTER  @Nwa_data; Preamble; Bytes;Binary( * ) 

297 

OUTPUT  @Fl;Preamblo;Byte8;Binary( *  ) 

296 

ASSIGN  @F1  TO  * 

299 

I 

300 

I  DISPLAY  REAL, IMAGINARY  PAIRS  ON  A  POLAR  PLOT 

301 

I 

302 

OUTPUT  @Nwa; "WAIT;CHAN2;WAIT;R£IP;WAIT; AUTO; " 

303 

OUTPUT  eNwa; "WAIT; CHANl ;WAIT;REIP; WAIT; AUTO; ENTO;” 

304 

J 

311 

1 

312 

1  SELECT  DELAY  VALUES  FOR  EITHER  OR  BOTH  CHANNELS 

313 

I 

317 

PRINT 

318 

INPUT  "INPUT  THE  ELECTRICAL  DELAY  IN  PICOSECONDS", Ed 

319 

OUTPUT  eNwa; "CHANl ; ELED" ; Ed; "pe; " 

320 

OUTPUT  eNwa ; " CHAN 2 ; ELED " ; Ed ; " pa ; ENTO ; " 

321 

1 

322 

1  ESTIMATE  THE  INITIAL  'N'  VALUE  FOR  COMPUTING  THE  CORRECT 

PHASE  ANGLE 

323 

4 

324 

OUTPUT  etlwa;  "CUAN2 ;  OELA;  WAIT;  " 

325 

OUTPUT  eNwa;"MARKr';Fi;";" 

326 

OUTPUT  eNwa; "OUTPMARK; " 

327 

ENTER  eNwa;Tau,0umitiy 

328 

N«INT(Tau*Fl ) 

329 

PRINT  "N,TAU  ",N,Tau 

330 

IF  N<0  THEN  N=0 

331 

IF  N>99  THEN  N-99 

332 

OUTPUT  eNwa; "WAIT;REIP;WAIT;AUTO; " 

334 

1 

335 

1  LOAD  THE  DATA  INTO  COMPLEX  ARRAYS  ON  THE  330  FOB  FURTHER  PROCESSING 

336 

I 

337 

OUTPUT  eNwa;  "WAIT;CHANl;FORM3;OUTPFORM; " 

330 

ENTER  eNwa  data; Preamble; Byte8;Sl 1 (• ) 

339 

OUTPUT  eNwa; "WAIT; CHAN 2 ; FORM3 ; OUTPFORM; " 

340 

ENTER  eNwa  data;Preaniblc;ByteB;S21(  *  ) 

341 

1 

342 

1  CALCULATE  THE  COMPLEX  DIELECTRIC  CONSTANT 

343 

1 

344 

BEEP  4000,. 2 

345 

C=3.E+10 

346 

RAO 

347 

PRINTER  IE  CRT 

348 

LINPUT  "DO  YOU  NEED  TO  SPECIFY  THE  NUMBER  OF  WAVELENGTHS  IN 

THE  SAMPLE?", A 

wlS 

349 

IP  AwlS-"N"  THEN  GOTO  351 

350 

INPUT  "INPUT  THAT  NUMBER", N 

351 

Nl6t-N 

353 

Anqlel'O. 

354 

lmrhB«0 
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356  INPUT  "TEST  VALUE  FOR  IMAG(RhE)  IS  O.C  DC  YOU  WANT  ANOTHER  VALUE?"  Alir.S 

357  IF  AiniS  =  "K"  THEN  GOTO  359 

358  INPUT  "INPUT  THE  CUTOFF  VALUE", Imrhs 

359 

360  FOR  1=1  TO  Np 

361  Rh6=(l  +  S21(I)*S21(I)-Sil(I)*Sll(I))/(2*S21(I)  ) 

362  IF  ABS(RhB)>l.E*6  77HEN  GOTO  382 

363  IF  .\BS( IMAG(RhB) )<Imrha  THEN  Rh6=CKPLX( REAL ( Rhs ) , 0 ) 

364  I,hs“ACS(Rh6) 

365  Lhs=CMPLX(REAL{Lh3) ,ABS(lKAG(Lhs) ) ) 

366  IF  IMAG(RhB)<0  THEN  LhB-CMPLX( 2*PI-REAL( Lhs ) , IMAG( Lhs ) ) 

367  Angle=REAL(Lhs) 

368  IF  I-l  THEN  GOTO  372 

369  IF  ABS(Ar.gle-Anglel)<4.5  THEN  GOTO  372 

371  N'Ntl 

372  LhB«=CMPLX(REAL(LhB)*N*2*PI,IMAG(Lh8)  ) 

373  1  PRINT  USING  "DDDD, 5( 2X,SDDD.DDD) , 2X,  DDD" ; I , Rhs , Lhe , Angle,N 

374  Anglel'Ajigle 

375  Rhs=Lh8*LhB 

376  Radfreq=2«PI*(Fi+(I-l)*(Ff-FiJ/(Np-l) ) 

377  Fact-C/ (Radf req*D) 

378  Epsp( I )=Fact"Fact*REAL(Rh8 ) 

379  Ep8pp(  I)=Fac"."Fact*IKAG(Rhs) 

380  1  PRINT  Ep3p( I ) ,EpBpp(I ) 

382  NEXT  I 

383  1 

384  :  PLOT  THE  RESULTS 

385  : 

386  BEEP  4000,. 2 

387  GINIT 

389  PRINTER  IS  CRT 

390  LINPUT  "INPUT  A  TITLE  FOR  THE  PLOT",TtlS 

391  PRINT 

392  PRINT  "MAXIMUM  VALUE  OF  EPSP=" ,MAX( Epsp( • ) ) 

393  PRINT 

394  PRINT  "MAXIMUM  VALUE  OP  EPSPP'" , MAX( Epspp( * ) ) 

395  PRINT 

396  INPUT  "INPUT'  THE  MINIMUM  AND  MAXIMUM  VALUES  OF  EPSILON  THAT  WILL  BE  PLOTTE 
0"  ,Waiin,Wmax 

397  INPUT  "INP'JT  THE  TICK  MARK  INCREMENT  IN  UNITS  OF  EPSILON" , Ytick 

390  CLEAR  SCREEN 

399  Xinax*=100"MAX(  1, RATIO) 

400  Yraax=100*MAX( 1, 1/RATIO) 

401  LONG  6 

402  CSIZE  1.7 

403  MOVE  Xmax*.3,Yniax 

404  LABEL  TtlS 

406  S0cllne5'"FILEr  "SNamefSS"  N»  "SVALS ( NlBt ) t "  DELAY*  "SVALS(Ed)6"  ps" 

408  MOVE  Xmax*.3,Ymax* .97 

409  lUlBEL  SecllneS 

411  DEG 

412  LUIR  90 

413  MOVE  Xniax*.04,Yinax*  .67 

414  LABEL  "EPS'" 

415  MOVE  XIMX*.04,Yra^oc•  .67 

416  LABEL  "EPS' '" 

417  LORG  4 

418  LDIR  0 

419  MOVE  Xmax* , 3 , . 55* Ymax 

420  LABEL  "FREQUENCY  (GHz)" 
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421  VIEWPORT  .  I’Xmax,  .  5‘XiEax,  .  8*Yrtax,  .95*Ymax 

422  WINDOW  7,  ll,Wroin,Wn'.ax 

423  AXES  1, Ytick,7,0, 1 , 1,2 

424  CLIP  OFF 

425  CSIZE  1.5, .5 

426  LORG  6 

427  FOR  I«e  TO  11 

426  MOVE  I,-.2»Ytick 

429  LABEL  USING  " « , K" ; ( 10* I ) / 1 . £♦ 9 

430  FOR  J-=2  TO  9 

431  Jtic)c«LGT(  J) 

432  MOVE  I-l+Jtick,Wroin 

433  PLOT  I-l+Jtick,  .04*  (Winax-Wnin) , -1 

434  NEXT  J 

436  NEXT  I 

437  LORG  e 

438  FOR  l=Wnln  TO  Winax  STEP  Ytlck 

439  MOVE  6.85,1 

440  LABEL  USING  " H ,SDDD . D" ; I 

441  NEXT  I 

442  PENUP 

443  FOR  I-l  TO  Np 

444  IF  ABS(EpBp(l) )<.001  THEN  GOTO  447 

445  X=>LGT{Pl*(I-l)*(Ff-Fl)/(No-l)) 

446  PLOT  X,Epsp(I) 

447  NEXT  I 

448  VIEWPORT  .  1  *Xjnax ,  .  5 *Xinax ,  .  6*  Ymax ,  .  75*  Yraax 

449  WINDOW  7,  ll,Wmlr.,Winax 

450  AXES  1, Ytlck, 7, 0, 1, 1,2 

451  CLIP  OFF 

452  CSIZE  1.5,  .5 

453  LORG  6 

454  FOR  1=8  TO  11 

455  MOVE  l,-.2*Ytick 

456  LABEL  USING  " # ,K" ; 10*1/ 1 . E+9 

457  FOB  J=2  TO  9 

458  JtlcK«LGT(J) 

459  MOVE  I-l«Jtlck,Winiri 

460  PLOT  l-l*Jtick,.04*(Wmax-Wmln),-l 

461  NEXT  J 

462  NEXT  I 

-.63  LORG  8 

464  FOB  I=Winln  TO  Wmax  STEP  Ytlck 

465  MOVE  6.85,1 

466  LABEL  USING  " # , SDDD . D" ; I 

467  NEXT  I 

468  PENUP 

469  FOR  I-l  TO  Np 

470  IF  ABS(EpBpp(I))>:.O01  THEN  GOTO  473 

471  X=LGT(Fl  +  (I-l)"(Ff-Fi)/(Np-l)  1 

472  PLOT  X,Epspp(I) 

473  NEXT  I 

474  KEY  LABELS  OFF 

475  PAUSE 

476  CLEAR  SCREEN 

47  7  KEY  LABELS  ON 

478  LOCAL  <?Nwa 

479  BEEP  4000, .2 

460  1  PRINT 

481  1  INPUT  "DO  YOU  WANT  TO  STORE  THIS  FILE  ON  A  BACKUP  TAPE?",An8S 
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482  !  IF  An8$="N"  THEN  GOTO  473 

483  I  copy  NaroefS  TO  Nanief$fi":  ,703" 

484  PRINT  "STOP  (1),  RETRIEVE  ANOTHER  FILE  (2),  SELECT  ANOTHER  DELAY  ( 3 V 

485  PRINT  "MAKE  ANOTHER  MEASUREMENT  (4)"  '  ' 

486  INPUT  sell 

487  IP  Sell-1  THEN  STOP 

488  IF  Sell-2  THEN  GOTO  222 

489  IF  Sell-3  THEN  GOTO  317 

490  OUTPUT  eNwa;"ENTO;" 

491  OUTPUT  eNwa;"WAIT;CHAN2;S21;CHANl;Sll;" 

492  Fd$-DATE$(TIMEDATE) 

493  FtS-TlME$(TIMEDATE) 

494  GOTO  263 

495  END 
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